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Abstract 

In recent years, the freely available Monte Carlo code REAS for modelling radio 
emission from cosmic ray air showers has evolved to include the full complexity of 
air shower physics. However, it turned out that in REAS2 and all other time-domain 
models which calculate the radio emission by superposing the radiation of the single 
air shower electrons and positrons, the calculation of the emission contributions 
was not fully consistent. In this article, we present a revised implementation in 
REAS3, which incorporates the missing radio emission due to the variation of the 
number of charged particles during the air shower evolution using an "end-point 
formalism" . With the inclusion of these emission contributions, the structure of the 
simulated radio pulses changes from unipolar to bipolar, and the azimuthal emission 
pattern becomes nearly symmetric. Remaining asymmetries can be explained by 
radio emission due to the variation of the net charge excess in air showers, which is 
automatically taken into account in the new implementation. REAS3 constitutes the 
first self-consistent time-domain implementation based on single particle emission 
taking the full complexity of air shower physics into account, and is freely available 
for all interested users. 
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1 Introduction 



High energy cosmic rays initiate extensive air showers when entering the 
Earth's atmosphere. The electromagnetic component of the air shower pro- 
duces radio emission which can be measured to study the characteristics of 
the cosmic rays. Radio detector arrays like LOPES p], [2] and CODALEMA 
[3] , [I] have verified that radio emission in the air is dominated by a geomag- 
netic effect, and they have studied correlations of the radio signal with shower 
parameters in great detail. With AERA [5] within the Pierre Auger Obser- 
vatory [6] and LOFAR [TJ, radio detection will be applied on larger scales 
and new "super-hybrid" techniques for measuring cosmic ray air showers will 
be employed. To understand the measurements and to learn more about the 
physics of cosmic rays using radio signals from air showers, a solid theoretical 
understanding of the radio emission process is needed. Presently, two major 
approaches exist, both of which are mainly based on geomagnetic effects [8]. 
On the one hand, there is a model based on the calculation of emission from 
the deflection of single particles in the geomagnetic field which was mainly 
investigated by Huege et al. [S], [TO], [H], [12], in particular with the imple- 
mentation in the REAS Monte Carlo code. The treatment of the radio emission 
with Monte Carlo techniques makes it straight-forward to couple it with de- 
tailed Monte Carlo air shower simulations. In case of REAS, this is done with 
CORSIKA [13] in which the needed particle distributions are histogrammed 
by COAST [12J. It is interesting to note that the original frequency-domain 
calculation of this implementation predicted spectra decaying to zero at small 
frequencies [9J, the time-domain implementation in the simulation code REAS, 
however, hitherto predicted spectra levelling off at a finite field strength for 
small frequencies, leading to essentially unipolar pulses [11], [12]. 
On the other hand, the macroscopic description of geomagnetic radiation 
(MGMR) developed by Scholten et al. [H] constitutes a modern implementa- 
tion of the approach for radio emission modelling by Kahn and Lerche in 1966 
|15j . Due to the Lorentz force, moving electrons and positrons are separated 
in the Earth's magnetic field. In MGMR, this is described by a time-variable 
net electric current in the electron-positron plasma, which is moving through 
the atmosphere with approximately the speed of light. The emission of an 
electromagnetic pulse is caused by these time-dependent transverse currents. 
Due to the dependence on the evolution of the charged particles in the shower 
with time, the MGMR model predicts a bipolar structure of the pulse shape 

m 

The differences in the results of both models were studied and led to the con- 
clusion that at least one model was not complete. The main reason for the 
different results of both approaches was that in REAS2 (and all other time- 
domain approaches based on single particle emission, e.g. [16], [UJ), emission 
due to the variation of the number of charged particles within the shower was 
not considered. A detailed discussion of this issue can be found in [18J to which 
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we kindly refer the reader for further details. To correct the implementation 
in REAS, these missing radiation contributions had to be taken into account. 
It should be stressed that the resulting implementation, based on the "end- 
point" formalism, consistently describes the radiation of the complete under- 
lying particle motion, not just "synchrotron" -like emission from the particle 
acceleration (see [19] for a discussion of the universal nature of the end-point 
formalism) . 



2 Simulation algorithm of REAS 



From REAS2 to REAS3, the simulation algorithm is unchanged and the air 
shower information is provided by CORSIKA and COAST in the same way 
as before. Therefore, we here only give a short overview of the technical im- 
plementation and for details, we kindly refer the interested reader to [T2] . 
First, the shower is simulated with CORSIKA using the air shower parameters 
of interest (such as primary energy, magnetic field, mass of primary, incoming 
direction, etc.). Using COAST, the information of the electrons and positrons 
is saved in histograms. These histograms contain information about the at- 
mospheric depth of the particle, the particle arrival time, the lateral distance 
of the particle from the shower axis, the particle energy and the particle mo- 
mentum direction. In the next step, REAS is generating individual electrons 
and positrons randomly according to the histogrammed distributions. These 
particles are then tracked analytically in the geomagnetic field. Note that the 
trajectories of the REAS simulation do not represent real trajectories of par- 
ticles, i.e. one long particle trajectory is represented by an ensemble of several 
shorter, unrelated trajectories in the code. The length of the trajectories is 
determined by a parameter A which is explained in chapter [5j 
In REAS3, the trajectory is distributed symmetrically around the position at 
which the particle is generated. 



3 Contributions due to charge variation 



In REAS2 [12], the radiation of single particles of an air shower is calculated 

as 



E(x, t) = e 



n-P 



^(l-P-nfR 2 



ret 



n x [{ft — (5) x 0\ 
{1-P-nfR 



ret 



where e indicates the particle charge, (3 = v(t) /c is given by the particle veloc- 
ity, R(t) = \R(t)\ describes the vector between particle and observer position, 
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n(t) = R(t)/R(t) is the line-of-sight direction between particle and observer, 
and 7 is the Lorentz factor of the particle. The index "ret" means that the 
equation needs to be evaluated in retarded time. 

The electrons and positrons emit radiation continuously along their track. To 
get a consistent description of radio emission in air showers, however, not only 
contributions due to the deflection of the particles in the magnetic field have 
to be taken into account, but also contributions due to the variation of the 
number of charged particles. 

REAS2 treats radiation processes only along the trajectories, but not at the 
end or the beginning of trajectories. Strictly speaking, this is equivalent to the 
situation that the particle arrives with the velocity v ~ c given by CORSIKA, 
enters the Earth's magnetic field where it is deflected and describes a short 
curved track and finally flies out of the influence of the magnetic field still with 
a velocity v ~ c. Figure [1] shows a sketch for such a particle trajectory. It is 




Fig. 1. Sketch of the the trajectory how 
it was implemented in REAS2. 




Fig. 2. Sketch of the the trajectory with 
start- and end-point as needed for a 
consistent description of radio emission 
from EAS. 



obvious that this is not describing the real situation in an air shower. Conse- 
quently, to complement the description in the Monte Carlo code, radiation at 
the beginning and at the end of each particle trajectory has to be calculated. 
If at a given atmospheric depth more particle trajectories start than end, i.e. 
the number of charged particles grows, this results in a net contribution. The 
same is true if the number of charged particles declines, i.e. more particle 
tracks end than start. Note that a net contribution occurs as well due to the 
change of the geometrical distribution caused by the spatial separation of the 
charged particles. 

In REAS3, a particle is treated as if it was created at rest and became "in- 
stantaneously" accelerated to v ~ c, flew on a short curved track through the 
Earth's magnetic field and got decelerated to rest again (see figure [2]). The 
acceleration process at the injection as well as the removal of an electron or 
positron takes place on time scales short in comparison with the frequencies 
of interest (V obsorvcd < 100 — 1000 MHz), i.e. 5t <C — - — . Hence, only the 

^observed 

time-averaged process is of interest, which will give a discrete contribution in 
contrast to the continuous emission along the curved particle trajectory. To 
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calculate contributions for the start-point, we consider the far field term of 
the general radiation equation and integrate over the injection time St. The 
"static" term of the radiation formula flTJ, the velocity field, can be neglected, 
because the "radiation" term completely dominates the signal for distances R 
relevant in practical applications. The relation for the retarded time used for 
transforming the integral from dt to dt' is derived from t' = t — R(t') / c. With 
dt = (1 — P ■ n)dt' the integral is solved as shown in the following calculation: 



E(x,t)dt 



c Jt 



n x [(n — (5) x (3} 



{1-/3- nfR 
*'i d I fix (fix p) 



dt 



CR Jt'r 



ret 



*i n x [(ft - P) x 13} 
(1-P-n) 2 



dt' 



cR Jt' dt 




1 - ■ ft) 

fix (fix p) 

n x (n x p) 



dt' 



cR Jt' 



Odt' 



[1-P-n) 

Likewise, one gains the electric field for the end-point of the trajectory: 

fix (fix p) 



(2) 



E(x,t)dt = - 



cR 



[1-P-n) 



(3) 



It is important to first transform the integration time into retarded time be- 
fore calculating the integral. The electric fields of the end-point^ZI and of the 
radiation along the track then have to be added to get a self-consistent imple- 
mentation for modelling radio emission in an air showers. 



4 Continuous vs. discrete calculation and incorporation of end- 
point contributions 



Adding the discrete endpoint contributions to the continuous contributions 
along the tracks may produce problems. The radiation associated with the 
end-points not only contains the emission due to the tangential acceleration. 
Due to the change in the direction of particle movement between the begin- 
ning and the end of the trajectory, radiation associated with the perpendicular 
acceleration, which was so far treated in the continuous description, is also 
contained. Combining the two descriptions therefore exhibits a risk of double- 

1 In the following we use the term "end-points" as a general term for both, start- 
points and end-points, as both are treated in the exact same way. 
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counting. In order to avoid such problems, it is preferable to change the cal- 
culation for the continuous contributions along the curved particle tracks to 
a discrete representation. Hence, the chosen representation in REAS3 is com- 
pletly discrete to ensure that the calculations are self-consistent. To describe 



v=0 




Fig. 3. Sketch of the trajectory with a discrete description and end-points. 

radio emission contributions along the particle trajectories in a discrete pic- 
ture, the trajectories of the particles are split in straight track fragments joined 
by "kinks" in which the velocity of the particles is changing instantaneously. 
A sketch of this description is given in figure [3j The instantaneous change of 
velocity at the kinks leads to radiation. With particle velocity 0i ■ c before and 
/?2 • c after the kink, radiation for one kink of the trajectory is: 



E(x,t)dt 




n x [(n — /3) x 0\ 



{1-13-nfR 

n x [n x (3 2 ) 

(1 - 02ft) 




(4) 



To verify that the continuous and the discrete calculations of emission con- 
tributions along the trajectories are equivalent and that both descriptions 
produce the same results, the REAS2 code was changed to calculate emission 
using the discrete approach of straight track segments connected with kinks 
(without the additional contributions of end-points). Analytically, the equiv- 
alence of the two approaches for the frequency domain has been shown by 
Konstantinov et al. [20]. Several tests with the REAS code have confirmed 
that the implementation of the discrete description is equivalent to the con- 
tinuous one as can be seen in figure HI The figure illustrates the equivalence 
for a vertical proton-induced air shower with primary energy of E p = 10 17 eV. 
The observer position is 100 m north of the shower core. The advantage of 
the discrete calculation is the consistency of the description of emission con- 
tributions along the tracks and emission at the endpoints which is making 
the incorporation of radiation at the endpoints canonical. To complement the 
former implementation in the Monte Carlo code with the emission due to the 
variation of the number of charged particles, it is therefore convenient to use 
the discrete description. In the discrete picture, contributions at the beginning 
and the end of a track are just kinks with @i = and f3 2 = 0, respectively. 
This self-consistent emission model has been incorporated in REAS3, taking 
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Fig. 4. Comparison of discrete and continuous description of radio emission (without 
endpoint contributions) for an observer 100 m north of shower core of a vertical 
10 17 eV air shower. The results are identical. 

the radiation at the beginning and the end of a trajectory as well as along 
the curved trajectory into account. The obtained results for the radio signal is 
discussed in section El but first the numerical stability will be demonstrated. 



5 Numerical stability 

As already mentioned in chapter [2} particle trajectories are represented by 
ensembles of several shorter, unrelated trajectories in the code. This ensures 
that the phase-space distribution of the analytically propagated particles stays 
consistent with the underlying particle distribution without having to treat 
energy losses during the propagation explicitly (cf. section 4 of [12] ). The 
length of the short segments is controlled by a parameter A. This parameter 
A determines the length of the short tracks into which the real trajectories 
of the particles are divided (it does not denote the sampling density of kinks 
on a track). If A is chosen inadequately large, the discrepancies between the 
particle distributions recreated in REAS and the distributions histogrammed 
in CORSIKA are getting too large. Thus, to avoid these discrepancies A has to 
be chosen small enough. In REAS2, it was recommended to set this parameter 
to 1 even though it is just a technical parameter, i.e. the result does not 
depend on the exact value of A as long as A is set small enough. Also, in REAS3 
the tracklength of the single short trajectories should have no influence on the 
physics results as long as it is chosen small enough. This cross check was 
made for a vertical proton-induced air shower with E p = 10 17 eV. As shown in 
figure [5] for an observer 100 m north of the shower core, the result converges 
with decreasing values of A. The same result is obtained for an observer 400 m 
north of the shower core which is displayed in the right column of figure |5j The 
comparison of both observer positions demonstrates that far away from the 
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shower core, the result converges faster than close to the shower core. For A = 
0.1 -^2 , a stable result is obtained in both cases. For observer positions close to 

cm* 7 1 




t [ns] t [ns] 

Fig. 5. Influence of different values of the technical parameter A. Left: observer 100 m 
north of shower core. Right: observer 400 m north of shower core. 

the shower core as for 100 m distance, it is also visible that for larger values of 
A, e.g. for A = 1.0 -^a , there is already a finite radio contribution for negative 
times. These unphysical contributions appear if partciles are created with 
symmetric trajectories around the place of particle generation and A is chosen 
inadequately large. For small enough values of A, however, the symmetric 
trajectories ensure a faster convergence of the result than the asymmetric 
trajectories used in RE AS 2. 

With smaller parameters of A, however, the high frequency noise increases 
as well as the computing time. This first effect is larger for observers far 
away as can be seen in figure EJ To optimize the calculation in REAS3, the 
path depth A of the electron and positron trajectories is therefore chosen as 
a dynamical parameter depending on the lateral distance of the observer. 
Therefore, it is recommended to set A to 0.1 at the core and increase this 
value linearly every 100 m as it is done by default in REAS3. The advantage of 
this implementation is on the one hand to gain stable results for all observers 
and on the other hand to avoid high frequency noise for observers at larger 
lateral distances. The difference in the recommended value of A in REAS2 and 
REAS3 is due to the fact that the emission model in REAS3 requires a precise 
description of the particle momenta during the propagation, which requires a 
more fine-grained treatment. 



6 Results 

In chapter \5\ we have verified that REAS3 is producing stable results and 
that the endpoints were implemented correctly. In chapter 16.11 we focus on 
the results obtained with REAS3 and in section 16.21 we discuss the influence 
of the charge excess in air showers on the radio emission. In both chapters, 
the radio emission was calculated for a proton-induced vertical air shower 
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with primary energy of E p = 10 17 eV. The shower itself was generated with 
CORSIKA 6.7 and COAST. The positions of the observers were chosen at sea 
level with different lateral distances and relative observer orientations to the 
shower core. For the simulations with REAS2 and REAS3, identical showers 
were taken. This is easily possible because the histogramming approach allows 
an easy seperation between the air shower modelling and the radio emission 
calculation. 



6.1 REAS2 vs. REAS3 



To study the changes introduced by the implementation of emission due to 
the variation of the number of charged particles, a simple, vertical shower 
geometry was chosen as specified above. The magnetic field was taken as hor- 
izontal with a field strength of 0.23 Gauss to get a geomagnetic angle of 90°. 
Figure |6] shows a comparison between unfiltered (i.e. unlimited bandwidth) 
pulses and frequency spectra of REAS2 and REAS3 for observers with dif- 
ferent positions with respect to the shower core. The spectra show the total 
field strength for two azimuthal observer directions: the thick line displays 
the total field strength for an observer east of the shower core and the thin 
line for an observer north of the shower core. It is clear that the strength of 
the pulses as well as the time structure of the pulses has changed, while the 
changes in the pulse amplitudes are dependent on the observer azimuth angle. 
The time structure of the pulses has changed from unipolar to bipolar. This 
can also be seen in the frequency spectra (lowest row of figure |6]) where in 
case of REAS3, the spectral field strengths drop to zero for frequency zero, 
as it was the case in the analytical implementation [9]. It can be argued from 
basic physical arguments that the spectral field strength has to drop to zero 
for small frequencies because the source of the coherent emission exists only 
for a finite time in a finite region of space and thus the the zero-frequency 
component of the emission, which corresponds to an infinite time-scale, can 
contain no power (cf. [2T]). 

While the emission pattern of REAS2 was azimuthally asymmetric, REAS3 
simulations are much more symmetric, as is expected for the given shower 
geometry. The increased symmetry can also be seen in the frequency spectra 
where the characteristics of the spectral field strengths with frequency are 
getting very similar for the two different azimuthal observer positions in case 
of REAS3. The remaining asymmetry will be discussed in section [6721 To get 
a more general impression of the changes in the amplitudes of the signal with 
REAS3, the lateral dependence of the emission is studied. Figure [7] shows this 
dependence for the unfiltered, full bandwith amplitudes. In REAS2, there was 
a much stronger signal for observers in the north or south than for observers 
in the east or west. In REAS3, the signal pattern is much more symmetric, 
but obervers in the eastern direction measure a higher absolute field strength 
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Fig. 6. Upper row: raw pulses for observers 100 m (left) and 400 m (right) north of 
shower core. Middle row: raw pulses for observers 100 m (left) and 400 m (right) east 
of shower core. Lowest row: frequency spectra for observers 100 m (left) and 400 m 
(right) east and north of shower core. For the raw pulses the east-west polarisation 
of the electric field is shown, whereas for the frequency spectra the total spectral 
field is shown. 



and observers in the west see a lower field strength than observers with other 
azimuthal positions. In general, the amplitudes of the field strength got lower 
from REAS2 to REAS3, while the changes for observers north and south of 
the shower core are much larger than for observers in the east or the west. 
To compare simulations with experimental data (e.g. of LOPES), the REAS 
simulations have to be filtered to a finite observing bandwidth. This is done by 
the helper application REASPlot which is included in the REAS3 package. In 
this paper, REASPlot was used with an idealised rectangular 43 MHz-76 MHz 
bandpass filter which can lead to acausal contributions at negative times due 
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Fig. 7. Lateral dependence with full bandwidth amplitudes. Left: REAS2. Right: 
REAS3. 

to the idealisation of the filter. In figure EJ the filtered pulses for observers 
with lateral distance of 100 m east and north of the shower core are shown. 
The quick oscillations are determined by the selected filter bandwidth. The 




Fig. 8. Filtered pulses (a simple rectangular 43-76 MHz bandpass filter is used) for 
observers with 100 m distance to the shower core. Left: the observer is north of the 
shower core. Right: the observer is east of the shower core. 

differences in the pulse strengths are much smaller for the filtered than for the 
raw pulses because the strongest changes between REAS2 and REAS3 occur 
at low frequencies (Please note that figure [6] is plotted on a log-log scale). In 
general, e.g. for different geometries, the differences in the amplitudes of the 
filtered pulses can be larger than in this example (cf. [22]). 
Again, the increased azimuthal symmetry of REAS3 is observable. To get an 
overall impression of the change from REAS2 to REAS3 and the influence 
of the net charge excess on the east-west asymmetry, we show the contour 
plots of the 60 MHz field strength in figure The ellipticity which is seen 
in REAS2 is replaced by a nearly azimuthally symmetric pattern in REAS3. 
The "clover" pattern seen earlier in the north-south polarisation component 
is no longer visible. In the vertical polarisation component (not shown here), 
there is no significant flux for either of the two simulations, as expected for 
a vertical shower. In the contour plots of REAS3, the east- west asymmetry 
is visible as well. The resulting REAS3 emission pattern can be interpreted 
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as a superposition of a circularly symmetric contribution with a v x B (in 
this case thus pure east-west) polarisation and a radially polarised emission 
contribution caused by the time- varying net charge excess. 
In summary, the incorporation of radiation due to the variation of the num- 
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Fig. 9. Contour plots of the 60 MHz field strength for emission from a 10 eV 
vertical air shower. From left to right: total field strength, north-south and east- west 
polarisation component. Contour levels are 0.1 ^Vm _1 MHz~ 1 apart. The closest 
position of the simulated observers to the shower core is 50 m. Upper row: REAS2. 
Lower row: REAS3 

ber of charged particles in the form of endpoint contributions results in a 
clear change from REAS2 to REAS3. The revised, self-consistent implemented 
model in the Monte Carlo code REAS3 predicts bipolar pulses with a mostly 
symmetrical emission footprint. In addition to these changes, a further emis- 
sion contribution arises from the variation of the net charge excess, which 
explains the remaining asymmetries. This contribution will be discussed in 
the following section. 

6.2 Discussion of charge excess emission 

The observed east-west asymmetry mentioned in section 16.11 arises from the 
fact that more electrons than positrons exist in an extensive air shower [23] . 
This net charge excess of order 10-20% leads to a contribution in the radio 
signal even in the absence of any magnetic field. Hence, in this section an air 
shower was simulated with CORSIKA for the geometry and primary charac- 
teristics as already mentioned at the beginning of this chapter, but in contrast 
to section 16. 1[ the magnetic field strength was set to Gauss. Although the 
showers calculated with B = 0.23 Gauss and B = 0Gauss do not represent the 
exact same particle distributions, we can interpret the predictions for the 
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shower with B = Gauss as the contribution due to the net charge excess of 
the shower used in section 16. 11 This allows us to test whether indeed the east- 
west asymmetry can be associated with the pure charge excess. In REAS2, 
radio emission is produced due to deflection of charged particles in the mag- 
netic field. Hence, there is no radiation without magnetic field. In contrast, 
REAS3 takes also emission due to the variation of the net charge excess into 
account. A radially polarised component for emission due to charge excess is 
expected, as seen also in the macroscopic approach [23]. The radiation pattern 
for a shower with B = Gauss is indeed radially polarised as illustrated in the 
contour plots of the 60 MHz field strength in figure [TUJ For the vertical com- 
ponent (not shown here) there is again no significant flux, as expected. Again, 
the closest observer position to the shower core is 50 m. Please note that the 
contour levels for the simulation without magnetic field are smaller than for 
the simulation with magnetic field, i.e., the relative field strength of the charge 
excess emission at 60 MHz is small in the distance range up to 200 m. To study 
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Fig. 10. Contour plots of the 60 MHz field strength for emission from a 10 17 eV 
vertical air shower without any magnetic field. Contour levels are 0.03 //Vm _1 MHz -1 
apart. The closest observer position to the shower core is 50 m. From left to right: 
total field strength, north-south and east-west polarisation component. The "spike" 
in the lower-left part of the contours is associated with noise in the simulation. 

the influence of the net charge excess emission on the overall radio signal, it 
is helpful to look at the polarisation vectors in the plane perpendicular to 
the shower axis. For the pure geomagnetic emission, the polarisation vectors 
at all observer positions point in the same direction as illustrated in the left 
sketch of figure (TTJ The right sketch illustrates the polarisation vector for the 



Fig. 11. Sketch of the polarisation vector. Left: uniform pattern as it is the case for 
pure geomagnetic emission. Right: radial pattern as it is the case for the net charge 
excess emission. 
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emission due to the variation of the net charge excess. The direction in which 
the vector points is changing with the observer position, following a radial 
pattern. Hence, for an observer in the east the total signal Se is given by 



Se — S, 



gin 



(5) 



where S gm is the pure geomagnetic contribution and S ce the net charge excess 
contribution to the signal. For an observer in the west, the total signal Sw is 
composed of 



S 



w 



Q _ Q 



(6) 



With the signals measured east and west from the shower axis, the signal for 
the pure geomagnetic emission and for the net charge excess can be calculated 
by 



gin 



-(Se + S 



w > 



and 



Sr. 



■(Se — Sw)- 



(7) 



To verify the assumption made in equations (EJ) and (jSJ), we calculated the 
signal of the charge excess as it is resulting from above for the shower with 
B = 0.23 Gauss and compared this with the emission of the shower with 
B = Gauss. Figure [121 illustrates that both pulses match. Therefore, the east- 
west asymmetry in the azimuthal emission pattern seen in the figures of section 





200 



Fig. 12. Comparison of emission due to net charge excess without magnetic field 
(solid line) and the calculated signal for the net charge excess (dashed line) from 
a shower with magnetic field. Displayed is in each case the east-west polarisation. 
Left: 100 m distance from shower core. Right: 400 m from shower core. 



16.11 is completely reducible to the emission of the net charge excess in an air 
shower. 

Finally, it is interesting to quantify the relative strength of the charge ex- 
cess emission with respect to the pure geomagnetic radio emission. Analyses 
studying the dependence of the radio signal on pure geomagnetic emission 
have already been done and have shown that there might be discrepancies 
between a pure geomagnetic model and the measured data [23]. The ratio of 
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the net charge excess signal and the pure geomagnetic radiation can be cal- 
culated from equation ((7j) to quantify the relative influence of the net charge 
excess. Figure [12] illustrates the ratio for the unfiltered full bandwidth ampli- 
tudes and the 43 to 76 MHz filtered bandwidth amplitudes for the east-west 
polarisation component. For the filtered case the ratio in the plot is shown 
only for observers up to 400 m lateral distance. The reason is that the fre- 
quency spectra for a vertical shower drop fast with increasing lateral distance 
as was also seen in figure |6j Consequently, in the used frequency range the 
signal is not anymore distinguishable from numerical noise at very large dis- 
tances. For the amplitudes of the unfiltered pulses the charge excess has more 
and more influence with larger distances, from a few % close to the core to 
around 90% of the pure geomagnetic emission at 1200 m. However, for the 
filtered pulses the ratio is almost constant over the whole range, at a level of 
~10%. It is important to clarify that the emission of the net charge excess 
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Fig. 13. Comparsion of the charge excess and the pure geomagnetic contribution on 
the radio signal of an vertical air shower for the east-west polarisation component. 
The lines may not represent the correct interpolation between the points due to 
the logarithmic scale of the x-axis. The ratio of the filtered data is only shown for 
distances up to 400 m, because then the signal is not anymore distinguishable from 
noise in this frequency range. 

occurs due to the variation of the number of charged particles and not due to 
Cherenkov-like emission. Both processes have been described in the pioneer- 
ing articles of Askaryan [26], [27], but the term "Askaryan radiation" is today 
generally interpreted as Cherenkov emission in dense media. For the inclusion 
of Cherenkov-like emission, a refractive index is needed which is not "unity" . 
REAS3 approximates the index of refraction to be unity so far. 



7 Conclusion 

Using an "end-point" formalism, we have implemented a fully consistent Monte 
Carlo model for radio emission from air showers, which in particular includes 
also the emission contributions due to the variation of the number of mov- 
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ing charges. For the implementation into the simulation code REAS, the 
description was changed from a continuous treatment to a discrete descrip- 
tion, because in the latter representation endpoint contributions can be added 
canonically. REAS3 produces stable results. The pulse shape has changed from 
unipolar pulses in REAS2 to bipolar pulses in REAS3. This is also apparent 
in the frequency spectra of REAS3 simulations which drop to zero for small 
frequencies. Several results illustrate the increased symmetry in the azimuthal 
emission pattern predicted by REAS3. Due to charge excess in extensive air 
showers it is evident that a small asymmetry has to remain. This effect was 
shown by calculating the pure geomagnetic and charge excess contributions 
from the asymmetries present in the full simulations and simulating radio 
emission in the complete absence of a magnetic field which agree well within 
numerical uncertainties. In the absence of a magnetic field, REAS3 predicts a 
radially polarized emission pattern. Due to the presence of this charge excess 
emission it is also obvious that the radio emission is not purely of geomagnetic 
origin and thus cannot be described by a pure v x B dependence, but that 
the signal polarisation depends on the exact observer position relative to the 
core. For a vertical shower, the relative strength of the emission due to the 
variation of the net charge excess with respect to the pure geomagnetic emis- 
sion increases from a few % close to the core up to ~90% at lateral distances 
of 1200 m for the unlimited bandwidth pulse amplitudes. This net charge ex- 
cess affects the east-west symmetry as well as the polarisation of the radio 
emission. For the filtered 43 MHz-76 MHz bandwidth amplitudes the relative 
strength of the net charge excess with respect to the pure geomagnetic emis- 
sion is constant at a level of approximately 10% for different lateral distances. 
Although the changes introduced with REAS3 are significant, many qualita- 
tive results obtained with earlier simulations are still valid (e.g., approximately 
exponential lateral distributions, a dependence of the lateral slope on X m3jX , 
the coherent scaling of the pulse amplitudes, and the fall-off of the frequency 
spectra to higher frequencies and larger distances). Certainly, the absolute 
amplitudes, pulse shapes and values of scaling parameters have changed with 
REAS3. With the implementation of emission due to the variation of the 
number of charged particles in REAS3, for the first time a self-consistent time- 
domain model exists which takes the full complexity of air shower physics as 
provided by CORSIKA simulations into account. The source code of REAS3 
will be freely available on request. REAS3 therefore constitutes a widely us- 
able radio simulation tool which can be employed to study radio emission from 
air showers initiated by arbitrary primary particles, using any of the avail- 
able hadronic interaction models, realistic atmospheric profiles and arbitrary 
shower geometries (including near-horizontal showers, since the atmosphere in 
REAS3 is treated with a curved geometry). 
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